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Abstract 

We have performed a holographic calculation of the hadronic contributions to the anomalous 

^ ; 

magnetic moment of the muon, using the gauge/gravity duality. As a gravity dual model of QCD 
' with three light flavors, we study a U(3)l x U(3)ij flavor gauge theory in the five dimensional 

^: 

AdS background with a hard-wall cutoff. The anomalous (electromagnetic) form factors for the 
, pseudo scalars, vr*^, rj and rj' , are obtained from the 5D Chern-Simons term of the gravity dual, 

' which correctly reproduce the asymptotic behavior of the form factor, dictated by QCD. We find 

o ' 

■ the total light-by- light contributions of pseudo scalars to the muon anomalous magnetic moment, 

. = 10.7 X lO^^'^, which is consistent with previous estimates, based on other approaches. 

a^ 
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> 
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Introduction 



One of most stringent tests of the standard model (SM) is provided by the measurement 
of the muon anomalous magnetic moment, a^, whose current precision is better than parts 
per million (ppm). Recent measurement of the {g — 2) value of the muon performed at 
the Brookhaven National Laboratory (BNL), 

= 11659208.0(5.4)(3.3) x 10"^° , (1) 

deviates by 2.2 — 2.7cr above the current SM estimate, based on e~^e~ hadronic cross sections. 
An improved muon {g—2) experiment has been proposed to achieve a precision of 0.1 ppm j^. 
The discrepancy between the SM estimate and the experimental value, if persists, might hint 
a new physics beyond the standard model. 

While the electroweak corrections can be calculated very precisely jsj, most uncertainties 
in the SM estimate of {g — 2) are coming from the hadronic corrections, which are essentially 
non-perturbative. The strong interaction contributions to the lepton magnetic moment 
consist of three pieces, the hadronic vacuum polarization, the higher-order hadronic vacuum- 
polarization effect, and the hadronic light-by-light (LBL) scattering. The contribution of the 
hadronic vacuum polarization is the leading correction of strong interactions to and 

has been recently calculated in lattice [4], but its present uncertainty is about 10 times larger 
than that of the current experiment. Fortunately, however, one can bypass the calculation of 
the hadronic- leading-order (HLO) correction js] and use the experimental result of e~^e~ 
(hadrons) (Fig.[T]), which is related to the hadronic vacuum polarization by the unitarity and 
analyticity of the diagrams. The higher-order hadronic vacuum polarization effect can be 
calculated quite accurately, once the hadronic vacuum polarization is obtained. It is found 
to be (5a^ = —101(6) x 10"^^ js]. Finally the hadronic light-by-light scattering correction 
(Fig. [2]) is the next-to-leading 0{a^) effect, but it is expected to be sizable to the current 
experimental accuracy, 6.3 x 10^^°. It is therefore absolutely needed to estimate its effect 
accurately to assess the SM deviation of the muon {g — 2). 

There have been several attempts to estimate the hadronic light-by-light scattering cor- 
rections, based on hadronic models or large Nc approximations 7, ^. For last two decades 
much improvement has been made in the LBL calculations despite a sign confusion joj. Cur- 
rently the positive sign is widel y q uoted for the pion-pole or pion-exchange contributions 



to the hadronic LBL correction 



IG 



ll| . which is dominant in the large N^. limit. Recently, 



however, a mistake related to the momentum conservation, which has not been considered 
seriously before, has been pointed out in the treatment of the 7r''7*7* form factor and 
is corrected properly jisl . 




Hadrons 



Quark loops 



FIG. 1: HLO correction (bottom) essentially involves a dressed propagator of the photon, which 
is related to the process (up) of annihilation into hadrons by unitarity and analyticity. 



The loop calculation of the LBL correction is logarithmically divergent in a hadronic 
model, where a constant vertex is used instead of full three-point form factors. This implies 
that the LBL correction is sensitive to the choice of the regulator or the cutoff A when 
we compare the model calculations with the data. In a model where the vector meson 
dominance (VMD) is adapted, however, the mass of vector meson naturally regulates the 
ultra-violet (UV) divergences. Among some models that regularize the hadronic uncertainty, 
the most popular ones are the meson dominance models like the VMD, the lowest meson 
dominance (LMD) and the LMD+V and so on, which fit among others the coefficient of 
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QCD axial anomaly for the low energy on-shell photon IJ, 

In this work we study the hadronic LBL contributions in holographic models of QCD, 
which naturally incorporate the vector meson dominance. Holographic models have been 



proposed recently for QCD 
string theory |l9l |. 
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18|, inspired by the gauge/gravity duality, found in the 



FIG. 2: Light-by-light correction (up) is supposed to be dominated by the diagrams, mediated by 
the pseudo scalar mesons (down). 



Several physical quantities of mesons and baryons such as their masses, couplings, and 



decay constants are calculated in holo grap 



agreement with the experimental data 
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lie models of QCD and found to be in a good 
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2l| . Encouraged by the success of holo- 



graphic QCD, we attempt to calculate the hadronic LBL contributions to the muon {g — 2). 
We consider in particular a hard-wall model of AdS/QCD, defined in a slice of five dimen- 
sional anti-de Sitter (AdS) spacetime, but our calculation can be easily applied to other 
holographic models of QCD. 



AdS/QCD and Light-by-Light scattering 

Solving QCD is very difficult, partly because at low energy hadrons, rather than quarks 
and gluons, are relevant degrees of freedom. Holographic QCD is to describe QCD directly 
with hadrons but in higher dimensions, prescribed by the gauge/gravity duality. The extra 
dimension is related to the energy scale of QCD, the boundary theory. In the gauge/gravity 
duality, the global symmetries associated with conserved currents of the boundary gauge the- 
ory become the gauge symmetries in the bulk, while the gauge-invariant boundary operators 
are mapped to bulk fields. 
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For a holographic model of QCD with three light flavors (up, down, and strange quarks) 
we consider a 5D f/(3)L x U{3)ji gauge theory in a slice of 5D AdS space-time, whose metric 
is given as, taking the AdS radius R= 1, 

1 



;{dx^dXfj_ — dz 



e < z < z„ 



(2) 



where the ultraviolet (UV) regulator e — >■ and an infra-red (IR) brane is introduced at z„ 



to implement the confinement of QCD 
by an action 



The model, known as AdS/QCD, is described 



S = j d^x^ Ti !^\DX\ 



+ Sy + S, 



CSi 



(3) 



where F is the field strength tensor of bulk gauge fields, A, dual of the QCD fiavor currents, 
and the 5D gauge coupling, g1 = 12tt'^/Nc {Nc is the number of colors). The bulk scalar 
X, which is bi-fundamental under the gauge group, is dual to the chiral-symmetry-breaking 
order parameter, (JlQr at the boundary, and its 5D mass, = A(A — 4), is related to 
the scaling dimension (A = 3) of the boundary dual operator. The covariant derivative 
is defined by DmX = OmX — iA^mX + iXA^iM^ where A^^r) denotes the U(3)L(ij) gauge 
fields. To correctly reproduce QCD flavor anomalies, one has to introduce a Chern-Simons 
term 221], given as 

^^^=2^/ N(Al)-cu5(A«)] , (4) 
where dujr^{A) = TrF^. Finally, for the anomalous U(1)a we introduce the bulk singlet Y, 
dual to G^^ (GG) of gluon flelds [23|, described by 

~1 



S- 



Y 



d^x^ 



,DY\' - -(r^'/det(X) + h.c.) 



(5) 



where the singlet Y has no mass term, since it is dual of dimension 4 operator, and k is a 
parameter to be flxed to give a correct mixing between rj and rj'. The bulk flelds A, X and 
Y, that we introduced, are related to towers of (axial) vector mesons [l8j and (pseudo) scalar 



mesons 2J], respectively, upon Kaluza-Klein reduction. There are also higher dimensional 
operators like or F^|Xp in the bulk action, but they are suppressed at low energy. 
(One could also introduce additional bulk flelds to describe higher-spin states, but they are 
irrelevant to our discussion.) 

According to the gauge/gravity duality, the classical bulk action becomes the generating 
functional for the one-particle irreducible (IPI) functions of the boundary gauge theory, 



once the source is identified as the bulk fields at the UV boundary. Therefore the hadronic 
LBL diagram, which is a four-point correlation function of the electromagnetic currents of 
quarks, can be easily calculated in holographic models of QCD, following the gauge/gravity 
correspondence. Since the quartic coupling for the bulk vectorial U(1)q^j^^ gauge fields, Vq^^ 
{Qem = |l + -^3), dual to the electromagnetic current of quarks, are absent in the bulk 
action ([3]), there are no IPI 4-point correlators for the electromagnetic currents. Therefore 
in AdS/QCD, where terms are suppressed, the LBL diagram is just given as a sum of 
IPI three-point functions, connected by intermediate states (FIG. |2]). The IPI three-point 
functions consist of three different types; the vector-vector-scalar correlators, the vector- 
vector-pseudo scalar correlators, and the vector-vector-axial vector correlators. Among them 
the vector-vector-scalar correlators are suppressed since they come from the higher order 
terms like while the vector-vector-pseudo scalar (axial vector) correlators are not 

suppressed as they are given by the bulk Chern-Simons (CS) term, Eq. (j4]). 

Once the LBL diagram is calculated, it is straightforward to evaluate the LBL corrections 
to the muon {g — 2). Our calculation is similar to that of Nyffeler [isl (see also 15|), where 
the pion contribution to the LBL correction was calculated most consistently. The only 
difference is that the full off-shell anomalous form-factors for pseudo scalars are derived 
from AdS/QCD rather than constructed to satisfy the QCD constraints. As we will see 
later, our form-factor derived from AdS/QCD does satisfy the asymptotic behavior for large 
and equal space-like photon momenta, obtained from perturbative QCD (pQCD), and also 
other asymptotic behavior derived from the operator-product expansion (OPE). 



vr*^, 7] and rj' form factor calculation in AdS/QCD 

We consider the anomalous pion form-factor first. Since the pions are decoupled from the 
strangeness fiavor, we need to consider up and down fiavors only, assuming equal mass. The 
vacuum solution of the bulk scalar field X is then written as {X) = j{mqZ + a z^)l2x2, where 
rriq and a correspond to the current quark mass and {qLlR), respectively. To analyze the 
correlation functions, we introduce the vector and axial- vector gauge fields V = {Al + Aji)/2 
and A = {Al — Ar)/2 and write X = {X) exp (227r"t"), where a = S denotes the abelian 
part, while a = a (= 1,2,3) are the SU(2) indices. (The generators = I2X2/2 and 
t°- = (J°-/2, a"' denoting Pauli matrices). Since the abelian component vr"^ is a part of rj, it 
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will be neglected for the discussion of the pion form-factor. 

In the axial gauge, V5 = = A^, the equations of motion for the transversal gauge fields 
are given as 



2 



ay 



0, 

0, 



(6) 
(7) 



where v{z) = rUgZ + az^ and V^(g, 2;) = V^{q)V{q, z) and A'^{q,z) = y4^(g)y4(g, z) are the 
4D Fourier-transform of vector and axial vector gauge fields, respectively. The normaliz- 
able solutions to Eq.'s (EI) and (I7j) will correspond to the vector and axial vector mesons, 
respectively, satisfying the boundary conditions V{q,e) = d^V^q, z)\z„, = (same for the 
axial vectors). The boundary conditions for the source fields (the non-normalizable modes) 
are on the other hand given as V{q, e) = 1 and dzV{q, z)\z^ = and same for A{q, z). The 
pion fields come from both vr" and the longitudinal components of the axial gauge fields, 
A\ = 9^0", which are related by equations of motion as, subjecting to boundary conditions 



r(g,e) = ar(g,^)i. 







7r"(g,e) 
dz ( -dzcl: 



-q'dz 



+ ^(7r'^-r) 
z-^ 

z^ 



0, 
0. 



(8) 
(9) 



We can easily read off the anomalous pion form-factor from the gravity dual, by taking 
the functional derivation of the bulk Chern-Simons action, Eq. (jl]) with the source field at 
the UV boundary, VQ^^(g,e). For arbitrary external photon momenta Q\ and (^2) the pion 
momentum being Qi ^ the form factor is found to be 

Nr. 



127r2 



i'iZm)JiQl,Zm)JiQ2 



dz dzip{z)J{Qi, z)J{Q2, z) 



(10) 



where ip{z), the wavefunction difference between 0" and vr", is the wavefunction of non- 
normalizable pion mode with boundary conditions, ^/'(e) = 1 and dz'4'{z)\z,n = ^"^^ J{Qi z) 
is the Wick-rotated expression of the solution V{q, z) defined below Eq. ([7]). The anomalous 
pion form-factor is shown in Fig. [3] for various kinematic regions. As argued in 25| , we find 
the pion form factor in AdS/QCD is in good agreement with the Brodsky-Lepage behavior, 
as shown in the Fig. [3l and also with the asymptotic behavior of pQCD for large space-like 
momenta of two photon legs, shown in the right panel of Fig. [31 We also note that our 



pion form-factor gives approximately the same magnetic susceptibility of quark condensate, 



obtained by OPE 
F . . 



26|. 
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FIG. 3: The anomalous pion form-factors 7*^7*7* (Q,0) from AdS/QCD (solid) and VMD model 
(dashed) are presented in the left panel. The right panel displays both FT^^*^*{Q^Qi) (lower two) 
and Ft^^*^*{Q,Q) (upper two) for both models. Notice that Ft^^*^*{Q,Q) grows as ~ asymp- 
totically. 



To calculate the anomalous form-factors for rj, rj', we need to consider the three-flavor 
case {Nf = 3) and introduce the term, needed for the U(1)a anomaly, included in the bulk 
action of flavor-singlet Y field, Sy, as r/ and t]' have mixing. The bulk gauge fields are now 
generalized to U(3)l x V{3)r and {X) = \{Mz + T.z^) with 



M 



( 



\ 







mg 
m, 



and S = o" X 1 



3x3 



(11) 



satisfying the similar equations as Eq.'s (jH]) and ([7]) except now v{z)/2 is replaced by the 
three-flavor vacuum solution (X). 

By solving the equations of motion for the pseudo scalars in the axial gauge, which are 
the phase fluctuations of bulk scalars, X and Y , and the longitudinal components of the 
axial gauge fields, the wave functions of r] and 77', that give the canonical kinetic terms 
for 7] and r/', can be found for an appropriate choice of {mri,K,) (FIG. S]) [23]. Once the 
correct wavefunctions are found, not only the amplitude of r], rj' 77 processes, but the 
whole momentum structure of the form factors can be surveyed through the same method 
used above to deduce the anomalous pion form-factor. Namely, the form factors are given 
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FIG. 4: The set of (m, k) which gives the canonical four dimensional kinetic term of the r/ and rj 
fields. Several masses are found for each k as a tower of the bound state of I^{J^^) = 
but only ground states are shown here. 

similarly to Eq. (ITOl) except that the non-normalizable pion wavefunction ip{z) is replaced by 

those non-normalizable wavefunctions of r/, r]' with different overall coefficients. As shown 

in FIG. ini their profiles are quite similar to those of pions. (Note also the anomalous form 

factors are correctly normalized to be consistent with the QCD anomalies.) Combining all 

these, we can now calculate the pseudo scalar contributions to the LBL corrections for the 

fl 

muon (^f — 2), following the recent work by Jegerlehner and Nyffeler [8|. 



Results and Discussion 

Numerical calculations in Table [T] has been performed with Mathematica 6. 0, using the 
adaptive MonteCarlo scheme. To utilize the method used in Ref. js], [3] for our calculation 
of meson exchange contributions, we have decomposed the source field of the vector gauge 
fields as follows [27]: 




(12) 
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FIG. 5: Pseudo scalar form factors Fpg^*^* (Q, 0) from the AdS/QCD are plotted. Q is a Euc 



momentum. Note that for Q = their offset values are well matched to the experiments 



28]. 



idean 



where ^pi^z) are the normalizable modes, corresponding to the (excited) rho mesons, with the 
boundary condition V'p(e) = and dz'^p{z)\z^ = and fp are their decay constants. Similarly 
we have expanded the bulk axial gauge fields A{q, z) for U(3) in terms of normalizable 
modes, which contain vr*^, r/, rj' and ai(1230), and towers of excited axial vector mesons. For 
our calculations we have set the parameters by Zm = 1/0.323, = 0.00222, = 0.04, 
a = 0.333^ and k = 35 where the energy scale is in GeV unit. We truncate for each photon- 
line the number of vector mesons in Eq. (1121) up to 4, 6, and 8 for our calculations. Beyond 
8th modes, the form factor changes very little, as shown in Fig. [6l We have shown our results 
in Table [H which are close to the recently revised value, a^^ = 9.9 x 10^^'^ in the LMD+V 
model by Nyffeler 13(], which correctly took into account of momentum conservation. 

To conclude, we have calculated the pseudo scalar contributions to the LBL corrections to 
the muon {g — 2) in an AdS/QCD model with three light flavors. Our holographic estimate 
gives results consistent with the recent estimate based LMD+V model. Our approach 
has a few parameters, which are highly constrained by low energy data of QCD, and is 
based on a principle, known as gauge/gravity duality, valid in the large Nc and large 't 
Hooft coupling (A) limit. Our result is therefore subject to l/N^ and 1/A corrections, which 
are at most 30%. 
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FIG. 6: Decomposing the source terms of vector gauge fields into the normalizable modes affects 
F^^,^, (Q,0) very little. The form-factor from VMD curve is dashed; yellow curve for 20 modes, 
blue for 60, and red for 150, respectively. 

TABLE I: Muon g - 2 results from the AdS/QCD. (10"^° factor should be multiphed to each 
number for a^.) 
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The holographic models of QCD generically show that the LBL diagram is given by a 
sum of IPI three-point functions, connected by intermediate states of pseudo scalars and 
axial vectors. For the LBL corrections, we have considered light pseudo scalars (vr^, 77, and 
Tj') only, since the axial vector meson, ai(1230), and its excited states are expected to be 
less important than light pseudo scalars. (A related work is in progress.) Finally we have 
calculated the full off-shell anomalous form-factors of light pseudo-scalars and found that 
they show the correct asymptotic behavior at large photon virtuality, given by pQCD. 
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